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We discuss the Taylor expansion approach to non-zero baryon chemical potential (/Xb) and present 
results on expanion coefficients of the pressure and energy density up to the 6-th order in j^b- 
Calculations have been performed with (2H-l)-flavor of improved staggerd fermions (p4fat3) and 
almost physical masses on lattices with temporal extent of four and six time slices. We use the 
expansion coefficients to construct the isentropic equation of state on lines of constant entropy per 
baryon number and various different hadronic fluctuations. Furthermore, we estimate the radius 
of convegence of our expansion, which can be seen as a method to analyze the structure of the 
QCD phase diagram. 
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Figure 1: Taylor coefficients of the pressure in term of the up-quark chemical potential. Results are obtained 
with the p4fat3 action on A^t = 4 (full) and Nx~6 (open symbols) lattices. We compare preliminary results of 
(2+l)-flavor a pion mass of niji k, 220 MeV to previous results of 2-flavor simulations with a corresponding 
pion mass of nip « 770 [^. 

1. Introduction 

A detailed and comprehensive understanding of the thermodynamics of quarks and gluons, 
e.g. of the equation of state is most desirable and of particular importance for the phenomenology 
of relativistic heavy ion collisions. Lattice regularized QCD simulations at non-zero temperatures 
have been shown to be a very successful tool in analyzing the non-perturbative features of the 
quark-gluon plasma. Driven by both, the exponential growth of the computational power of re- 
cent super-computer as well as by drastic algorithmic improvements one is now able to simulate 
dynamical quarks and gluons on fine lattices with almost physical masses. 

At non-zero chemical potential, lattice QCD is harmed by the "sign-problem", which makes 
direct lattice calculations with standard Monte Carlo techniques at non-zero density practically 
impossible. However, for small values of the chemical potential, some methods have been success- 
fully used to extract information on the dependence of thermodynamic quantities on the chemical 
potential. For a recent overview see, e.g. 

2. The Taylor expansion method 

We closely follow here the approach and notation used in Ref. |^]. We start with a Taylor 
expansion for the pressure in terms of the quark chemical potentials 




(2.1) 



The expansion coefficients c"'j'^{T) are computed on the lattice at zero chemical potential, using 
stochastic estimators. Some details on the computation are given in [|3|, Details on our cur- 
rent data set and the number of random vectors used for the stochastic random noise method are 
summarized in Table |l[ 

In Fig. d we show results on the diagonal expansion coefficients with respect to the up-quark 
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Table 1: Details on the calculation: The columns give from left to right the temperature values, the number 
of evaluated configurations, the number of trajectories by which these configurations are separated and the 
number of random vectors used for the evaluation of the traces, for A^t = 4 and 6, respectively. 

chemical potential up to the six order (c"'q q with n = 2,4,6). They have been evaluated on con- 
figurations generated for a calculation of the equation of state obtained with (2+l)-flavor of 
the p4fat3 action [^] and a pion mass of mj^ « 220 MeV. Here the full symbols are from = 4 
lattices, while the open symbols denote results from = 6 lattices. The spatial lattice size has 
been chosen to be = AN^, which is known to be large enough in order to suppress finite volume 
effects. Considering cut-off effects, we find that they are small and of similar magnitude as those 
found for the trace anomaly [Q]. This was already anticipated by the analysis of the cut-off correc- 
tions that arises in the free gas limit [Q]. Similar results for the expansion coefficients have been 
also obtained with the asqtad ^ and the standard action |^. 

We also compare our preliminary results for (2-1-1 )-flavor QCD and a pion mass of niji 
220 MeV, with previously obtained result of 2-flavor QCD and irijt « 770 MeV [|] (also p4fat3 
and Nt = 4). It is apparent from Fig. [l] that the critical temperature for these two particular sets 
of lattice parameter differ substantially. In fact, the transition temperature decreases from about 
225 MeV for the heavier mass calculations to about 200 MeV for the lighter mass calculations. 
Note, that those Tc values are the = 4 values, which of course are still influenced by the finite 
lattice spacing. Furthermore, we find from Fig. |l| that the quark number fluctuations of second, 
fourth and sixth order, which are related to these expansion coefficients, increase with decreasing 
quark mass. 

Alternatively to the quark chemical potentials one can introduce chemical potentials for the 
conserved quantities baryon number B, electric charge Q and strangeness S {piB,Q,s)^ which are 
related to ^u4,s via 

12 11 11 

llu = ^^iB + ^llQ, lid = -^iB - ::llQ, IX, = -IXb--I^q-Hs- (2.2) 
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Figure 2: Taylor coefficients of the pressure, the energy and entropy density with respect to the baryon 
chemical potential. Results are obtained with the p4fat3 action on A^^ = 4 (full) and Nr = 6 (open symbols) 
lattices. 



BOS 

By means of these relations the coefficients 'jj^ of the pressure expansion in terms of IIb,q,s are 
easily obtained, in analogy to Eq. 11 



(2.3) 



i,i,k 



For the rest of this article we will restrict ourselves to the case of Hq = Hs = 0, thus we will 
suppress in the following the indices that are related to these chemical potentials. From the pressure 
we immediately obtain the baryon number density n^, which is given by the derivative of p/T^ 
with respect to the baryon chemical potential /ig and can be expressed in term of the expansion 
coefficients , we have 

(2.4) 

n=2 



Using standard thermodynamic relations we can also calculate the expansion coefficients of the 
trace anomaly or equivalently the difference between energy density and three times the pressure, 



E-3p 

f4 



n=0 



(2.5) 



where the expansion coefficients are given by 



dT 



(2.6) 



Combining Eqs. 2.3, 2.5, and 2.6 we then obtain the Taylor expansions for the energy and entropy 
densities [fi^ 



£(3cf(r)+c-f(r)) 
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n=0 

j_ ^ e+p- jJ-BnB 



n=0 
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(2.7) 
(2.8) 



4 



Non-zero density QCD by the Taylor expansion method 



Chuan Miao and Christian Schmidt 



At present, we calculate the expansion coefficients from the coefficients , in accordance with 
Eq. by performing the T derivative numerically, which introduces a small systematic error. 

In Fig. § we show the second, fourth and sixth order expansion coefficients of the pressure, 
energy density and entropy density as given in Eqs. and |2.8| , obtained with the p4fat3 action. 
Full symbols are from Ni=A lattices, while the open symbols denote results from N^ = 6 lattices. 
We again find small cut-off effects, however, higher order derivatives of pressure, energy density 
and entropy density with respect to /ig are still very preliminary, as the error bars are large. This 
is especially true for the results from = 6 lattices. Nevertheless, the overall pattern of the 
coefficients is in agreement with expectations based on an analysis of the singular behavior of the 
free energy, making use of an appropriate scaling Ansatz. 

We find that the magnitude of the coefficients is decreasing drastically with increasing order, 
for all analyzed temperatures. Thus an approximation of the equation of state for small baryon 
chemical potential by means of a fourth or sixth order expansion seems to be justified. In general, 
an analysis of the radius of convergence of such a Taylor series is of great interest for an analysis 
of the QCD phase diagram, since the radius of convergence is bounded by the location of the QCD 
critical point as well as by any first order phase transition line (see Sec. 



3. The isentropic equation of state 



By using the Taylor expansion coefficients of the baryon number (Eq. 2.4) and entropy density 



(Eq. 2.8), we can compute the ratio of entropy per baryon number as function of T and jUg. The 0-th 
order coefficient for the entropy density has been taken from [^. Solving numerically for a constant 
ratio of entropy per baryon number, s/ub, we determine isentropic trajectories in the (r, /Xgj-plane. 
These trajectories are relevant for the description of matter created in relativistic heavy ion colli- 
sions. After equilibration the dense medium created in such a collision will expand along lines of 
constant entropy per baryon. It then is of interest to calculate thermodynamic quantities along such 
isentropic lines. 

We find that isentropic expansion at high temperature is well represented by lines of constant 
Hb/T down to temperatures close to the transition, T ~ 1.2ro. In the low temperature regime we 
observe a bending of the isentropic lines in accordance with the expected asymptotic low temper- 
ature behavior. The isentropic expansion hnes for matter created at SPS correspond to ~ 45 
while the isentropes at RHIC correspond to s/ns ~ 300. The energy range of the AGS which also 
corresponds to an energy range relevant for future experiments at FAIR/Darmstadt is well described 
by s/ng ^ 30. These lines are shown in Fig. ^ (left) together with data points characterizing the 
chemical freeze-out of hadrons measured at AGS, SPS and RHIC energies. These data points have 
been obtained by comparing experimental results for yields of various hadron species with hadron 
abundances in a resonance gas [ [IT| , [I^ ]. The solid curve shows a phenomenological parameteri- 



zation of these freeze-out data [12]. In general our findings for lines of constant s/ub are in good 
agreement with phenomenological model calculations that are based on combinations of ideal gas 
and resonance gas equations of state at high and low temperature, respectively [|l^, 14]. 



Results shown in Fig. ^ are based on a fourth order expansion of the pressure. We find, 
however, that the truncation error is small, i.e. the results change only little when we consider also 
the sixth order term in /ig. In accordance with the good convergence of our results, we find, that all 
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Figure 3: On the left panel we show Isentropic trajectories in the {T — /iB)-diagram, corresponding to 
s/nB ~ 300,45,30 respectively. Open symbols are form A^^ = 6, obtained by a 4th order Taylor expansion 
of the pressure. Filled symbols are from A^^ = 4 calculations. We also show freeze-out data, as well as a pa- 
rameterization of the freeze-out curve from [p^. Corresponding free gas limits are IIb/T = 0.21, 1.40,2.13 
respectively and are indicated by solid lines. On the right panel we plot the the ratio of pressure and energy 
density along those trajectories. 



trajectories shown in Fig. || (left) are well within the radius of convergence of the Taylor series. At 
present we estimate the radius of convergence of the pressure series to {hb/T)™^>2.1. The cut-off 
effects can be estimated by comparing open and full symbols. 

We now proceed and calculate energy density and pressure on lines of constant entropy per 
baryon number using our Taylor expansion results up to ^(/i^). Again we take the 0-th order 
coefficients from We find that both quantities obtain corrections of about 10% at AGS (FAIR) 
energies (s/ns = 30) and high temperatures. The dependence of s and p on s/ub cancels to a large 
extent in the ratio p/s, which is most relevant for the analysis of the hydrodynamic expansion of 
dense matter. This may be seen by considering the leading ^{jJ-g) correction, 



P _ 1 1 £o - 3po 
£ 3 3 £o 



C2 £2 



£() — ^Pq £() 



(3.1) 



In Fig. ^ (right) we show p/e as function energy density along our three isentropic trajectories. 
The softest point of the equation of state is found to be {p/e)mm — 0.07 — 0.09, for Nt = 4 and 
6 respectively. Within our current numerical accuracy it is independent of s/ub- Similar results 
for the asqtad action have been obtained in However, as our data is preliminary, the analysis 
clearly suffers from insufficient statistics, which is in particular true for our = 6 results. 



4. Hadronic Fluctuations 



Quark number fluctuations are related to the derivatives of QCD partition function with respect 
to the quark chemical potentials by the fluctuation-dissipation theorem. The Taylor expansion co- 
efficients cl'j'^, as defined in Eq. 2.1, can thus be directly interpreted as quark number fluctuations 
at /I = 0. However, quark fluctuations can not be detected directly in experiments due to confine- 
ment. Therefore we will consider fluctuations in terms of hadronic quantum numbers, i.e. baryon 
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Figure 4: Quadratic and quartic fluctuations of baryon number, electric charge and strangeness, normalized 
by their corresponding Stefan-Boltzmann value. The results on A^^ = 4 lattices (open symbols) and Nt: — 6 
lattices (full symbols) are in good agreements. 



number B, electric charge Q and strangeness S, which are more easily obtained by experiment. 

BOS 

These fluctuations are related with the hardonic Taylor expansion coefficients c- 'J^ , as given in 
Eq. U. 

In general, the quadratic fluctuations X2 chemical potentials can be obtained from the 



second order coefficient cf 



2VT^d{^x/Ty 



I^B.Q.S=0 



(4.1) 



where 5N = N — (N) denotes the normalized net-density and (. . .)o indicates that the expectation 
value has been taken at IJ.b,q,s = 0- Under such conditions, baryon number, electric charge and 
strangeness vanish, and we have 5N = N. We define 



1 



:2cf 



Similarly, we obtain the quartic charge fluctuations by 

1 



{{N'x)o-3{Ni)l) 



24cf , 



and correlations among two conserved charges by 

1 



All 



X22 



1 



{{NxNy)o-{Nx)o{Ny)o)=c'1{ , 
{{NlN^)o - {N^)o{N})o - 2{NxNy)1) = 4cf, 



Y 
22 ' 



(4.2) 

(4.3) 

(4.4) 
(4.5) 



where X,y e {B,Q,S}. 

In Fig. ^ we show results for quadratic and quartic fluctuations of B, Q and S. The quadratic 

BOS BOS 

fluctuations X2 rise rapidly in the transition region where the quartic fluctuations X4 show a 
peak. The peak height is more pronounced for the baryon number fluctuations than for fluctuations 
of the strange quarks. 
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We compare the results obtained on lattices with temporal extent A^^ = 4 and 6. We notice that 
they are in general compatible with each other, especially in the high temperature phase, where both 
quadratic and quartic fluctuations approach the Stephan-Boltzmann limit quickly. The transition 
temperature has been previously determined to be = 202 MeV and 196 MeV on A^, = 4 and 
6 lattices respectively [p3|]. We thus conclude that at temperatures of about LST^ and higher, 
quadratic and quartic fluctuations of B, Q and S are well described by the ideal massless quark gas. 

At low temperature, hadrons are the relevant degrees of freedom. The hadron resonance gas 
(HRG) model has been shown to provide a good description of thermal conditions at freeze-out. We 
thus compare the fluctuations in the low temperature phase with a HRG model, where we include 
all mesons and baryons with masses smaller than 2.5 GeV from the particle data book. 

In Fig. ^ we show the ratio of quartic and quadratic fluctuations for B, S and Q. In the 
HRG model, Xa /X2 easily obtained in the Boltzmann approximation, which is valid for a dilute 
baryonic gas in the temperature range of interest. One finds that all details on the hadron mass 
spectrum and temperature dependence cancel and the result is a constant, given by the unit square 
of the baryonic charge (one for all baryons). This is in fact reproduced by the lattice results shown 
in Fig. H (top left). 

The ratio of quartic and quadratic fluctuations for S and Q are more complicated even in the 
Boltzmann limit, since hadrons with different electric/strange charge give rise to different contribu- 
tions to the corresponding fluctuations. For strangeness fluctuations, shown in Fig. ^ (bottom), the 
Boltzmann limit is still a good approximation; but for electric charge fluctuations, the pion mass 
plays an important role. In order to check for the sensitivity of electric charge fluctuations on the 
pion mass, we show in Fig. ^ (top right) results of a HRG model calculation with physical pion 
masses and without the pion sector, i.e. for infinitely heavy pions. 

In Fig. ^, we show the various correlations Xii^' Xn Xn normalized to quadratic fluc- 
tuations Xi and X2 respectively. The results from A^^ = 4 and 6 lattices agree with each other 
very well, and they are compared with the HRG model in the low temperature phase and Stephan- 
Boltzmann limit in the high temperature phase. We find that the correlations reproduce the quali- 
tative behaviour of the HRG model below Tc and again start to agree with the free gas predictions 

forr>i.5r<. 



5. Radius of convergence 

Whenever an observable at jLt^ > is approximated by a truncated Taylor series, expanded 
around /ig = 0, an analysis of the convergence behavior of the series is, of course, mandatory. In 



addition, the convergence radius of such a series, as e.g. given in Eq. for the pressure, is of 
great interest for the structure of the QCD phase diagram. By definition, the radius of convergence 
is the distance to the closest singularity in the (r,/XB)-plane. It is bounded by the QCD critical 
point as well as by any line of first order transitions and could thus provide an interesting method 
to determine the QCD critical point. This method was proposed in Ref. [^] and was applied in the 
case of standard staggered fermions in Ref. 

One way to evaluate the radius of convergence p is by means of the formula 



p(r) = lim p„(r), with Pn{T) = Jcm/cl^^{T). (5.1) 
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Figure 5: Ratio of quartic and quadratic fluctuations of baryon number (B), strangeness {S) and electric 
charge {Q). The two curves of the HRG model in the top right figure correspond to charge fluctuations with 
physical pions (upper) and infinitely heavy pions (lower curve). 



In practice, we can not perform the limit « — > oo, instead we consider the first few approximations 
Pn{T), with «<6. The second approximation for p(T) is given by P2{T) = y^'cf (r)/cf (T) = 

VT^/y^^f/^- We see that X4/X2 shown in Fig. ^ is closely related to the radius of con- 
vergence. If this quantity develops a peak, which increases above the resonance gas value of 
X4 /Xi ~ 1' might expect the critical point to be found at (/iB/r)™'<\/l2. Of course, higher 
order approximations for the radius of convergence are needed to substantiate such an estimate and 
to establish the existence of a critical point. 

Our results are at present not conclusive, as the = 6 data shown in Fig. ||(left) suffers from 
insufficient statistics and temperature resolution. Note, however, that the peak in the baryon number 
fluctuation ratio xf IXi will develop in the chiral limit as can be seen from a scaling analysis of the 
free energy and has recently also shown in a chiral model calculation [[l^]. 

6. Conclusions 

We have calculated corrections to the equation of state arising from a non-zero baryon chemi- 
cal potential, by means of a Taylor expansion of the pressure. Within this framework we calculated 
the isentropic equation of state along lines of constant entropy per baryon number (s/ub) for RHIC, 
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Figure 6: Pairwise correlatios of the conserved charges baryon number (B), electric charge (Q) and 
strangeness (S) as function of the temperature, normaHzed to the quadratic fluctuations of B and Q re- 
spectively. 



SPS and AGS (FAIR) energies. Within our current, preliminary, analysis we find the softest point 
of the equation of state to be independent of s/ns- 

Furthermore, we have analyzed the quadratic and quartic fluctuations of baryon number, elec- 
tric charge and strangeness, as well their ratios. We find these quantities to be in good agreement 
with the free gas results at temperatures of T > 1.57^. Below T^, qualitative features of the res- 
onance gas are reproduced. This ratio for the baryon number is closely related to the second 
approximation of the convergence radius of the Taylor series of the pressure with respect to the 
baryon chemical potential. 
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